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We give short introduction to the Analytic Perturbation Theory (APT) [T] in QCD, 
describe its problems and suggest as a tool for their resolution the Fractional APT 
(FAPT) [2l |3]. We also describe shortly how to treat heavy-quark thresholds in 
FAPT and then show how to resum perturbative series in both the one-loop APT 
and FAPT. As applications of this approach we consider the Higgs boson decay 
H° — > bb, the Adler function D{Q^) and the ratio R{s) in the Nf = 4 region. Our 
conclusion is that there is no need to calculate higher-order coefficients d„>5 if we 
are interested in the accuracy of the order of 1%. 

1 Basics of APT in QCD 

In the standard QCD Perturbation Theory (FT) we know that the Renormalization 
Group (RG) equation das[L]/dL ~ ~al~ . . . for the effective coupling as{Q^) — as[L]//3/ 
with L = ln(gVA2), = 6o(iV/)/(47r) = (11 - 2Nf/S)/{A7rj^ Then the one-loop 
solution generates Landau pole singularity, as[L] = l/L. 

In the Analytic Perturbation Theory (APT) we have different effective couplings in 
Minkowskian (Radyushkin [3], and Krasnikov and Pivovarov f5]) and Euclidean (Shirkov 
and Solovtsov [l]) regions. In Euclidean domain, -q^ ^ ^ L ^ InQV^^, APT 
generates the following set of images for the effective coupling and its n-th powers, 
{A,i[i]}„gpj, whereas in Minkowskian domain, = s, Lg = Ins/A^, it generates another 

*This work is partially done in collaboration with S. Mikhailov (JINR) and N. Stefanis (ITP-II RUB). 
^We use notations f{Q^) and f[L] in order to specify what arguments we mean — squared momentum 
or its logarithm L = ln(Q^/A^), that is f[L] = f{Is? -e^) and is usually referred to N j = 3 region. 
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set, {2l„[Ls]}„gN- APT is based on the RG and causality that guaranties standard per- 
turbative UV asymptotics and spectral properties. Power series c?ma^[L] transforms 
into non-power series X^m^™-^™!^] APT. 

By the analytization in APT for an observable f{Q^) we mean the "Kallen-Lehman" 
representation 

Then in the one- loop approximation (note pole remover (e^ — in (|2ap ) 



a + g2 L - 1 



ML] = / T^'^^ = T-T^> (2a) 



oo 



2ti[L,] = / ^ rfa = - arccos , (2b) 

7^ V + 



whereas analytic images of the higher powers (n > 2, n e N) are: 



%n[Ls\J {n-l)\ V dL) V2li[i 



(3) 



In the standard QCD PT we have also: 

(i) the factorization procedure in QCD that gives rise to the appearance of logarithmic 
factors of the type: [L] L; El 

(ii) the RG evolution that generates evolution factors of the type: B{Q^) = [Q'^) / Z [ji^)]^ 
B{iJ,^), which reduce in the one-loop approximation to Z{Q'^) ^ [L] with i/ = 7o/(26o) 
being a fractional number. 

All these means we need to construct analytic images of new functions: , L™, .... 

In the one-loop approximation using recursive relation ([2]) we can obtain explicit 
expressions for Au[L] and 2li/[i]: 



1 F{e-\l-i.) „ _ sin [(^ - 1) arccos {L/Vt:^ + L^)] 

Here i^(z, v) is reduced Lerch transcendental function, which is an analytic function in 
ly. Interesting to note that Ai,[L\ appears to be an entire function in i/, whereas 21. [L] 



^ First indication that a special "analytization" procedure is needed to handle these logarithmic terms 
appeared in [6], where it has been suggested that one should demand the analyticity of the partonic 
amplitude as a whole. 
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is determined completely in terms of elementary functions. They have very interesting 
properties, which we discussed extensively in our previous papers [SjrT. Here we only 
display graphics of Ai, [L] and 21^ [L] in Fig. [TJ one can see here a kind of distorting mirror 
on both panels. 
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Figure 1: Graphics of A,y[L] (left panel) and 2ty[L] (right panel) for fractional v [2, 3]. 

Construction of FAPT with fixed number of quark flavors, Nf, is a two-step procedure: 
we start with the perturbative result [as((5^)] ^ , generate the spectral density Pv{cr) using 
Eq. and then obtain analytic couplings A,^[L] and 2t,y[L] via Eqs. Here Nf is fixed 
and factorized out. We can proceed in the same manner for A'^^-dependent quantities: 
[a,iQ^;Nf)Y ^ p.(a;7V/) = p.[L.;iV/] = pA<y)/^f ^ A,[L;Nf] and 2l,[L;iV/] - 
here Nf is fixed, but not factorized out. 

Global version of FAPT, which takes into account heavy-quark thresholds, is con- 
structed along the same lines but starting from global perturbative coupling [a|'°''((5^)] " , 
being a continuous function of due to choosing different values of QCD scales A/, 
corresponding to different values of Nf. We illustrate here the case of only one heavy- 
quark threshold at s = m|, corresponding to the transition A^^ = 3 — > Nf — 4. Then we 
obtain the discontinuous spectral density 

Pf"!'^) = PtMi^] ^e{L„< U) p„ [L,; 3] +e{U< L,) p„ [L, + A4; 4] , (5) 

with La = ln(cr/A|), Lf = ln(^m^/A|j and A/ = In (^A|/A^^ for / = 4, which is 
expressed in terms of fixed-flavor spectral densities with 3 and 4 flavors, p„[L;3] and 
pn[L + A4;4]. However it generates the continuous Minkowskian coupling 

2trii.] = (?(£.<i4)(2l,[£.;3]-2l,[L4;3]-f 2l,[L4 + A4;4]) 

+ e{U<L,)%^[Ls + \i-A]- (6) 
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and the analytic Euclidean coupling (for more detail see in [7]) 



p^[L^;3] - p^[L„ + A4;4] 



1 + e 



L-L„ 



dL„ 



(7) 



2 Resummation in the one-loop APT and FAPT 

We consider now the perturbative expansion of a typical physical quantity, like the Adler 
function and the ratio i?, in the one-loop APT 



V[L\ 

n[L] 



= dn 



'An[L] 



(8) 



We suggest that there exist the generating function P{t) for coefficients dn — dn/di 



P{t)e-^dt with 



P(t) = 1 



(9) 



To shorten our formulae, we use the following notation ((/(t)))p(t) = J^f{t)P(t)dt. 
Then coefficients dn = di {{t"~^P{t) and as has been shown in [8J we have the exact 
result for the sum in IH]) 



V[L]\ ^ 

nw 



da + di 



{{Ai[L~t]))p^t) 
((2li[L-t]))p(t) 



(10) 



The integral in variable t here has a rigorous meaning, ensured by the finiteness of the 
couplings Ai[t] < 1 and 2li[t] < 1 and fast fall-off of the generating function P{t). 

In our previous publications [71 [3] we have constructed generalizations of these results, 
first, to the case of the global APT, when heavy-quark thresholds are taken into account. 
Then one starts with the series of the type (O, where An[L] or 2l„[L] are substituted 
by their global analogs, An°^[L] or 21^'°'' [L] (note that due to different normalizations of 
global couplings, ~ An[L]//3f, the coefficients dn should be also changed). The 

most simple generalization of the summation result appears in Minkowski domain: 



7^«""'[L] do + di{{e{L<U) A42ti[i]+2t] 
+ di{{e{L>U)^ 



P3 ■ 



Pit) 



L + A4-— ;4 

P4 ■ 



» 



Pit) 



(11) 
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where A4^^[t] L4 + X4- i/A; 4 - 2t^ L3 - t/Ps; 3 . 

The second generalization has been obtained for the case of the global FAPT. Then 
the starting point is the series of the type X^i^o '^^ ^n+v I^l ^^'^ result of summation 
is a complete analog of Eq. ([TTjl with substitutions 



P{t) ^P,{i)= P 



1 



1 



(12) 



do =^ do2l^[i], 2ti[L-t] => %i+u[L-t\, and A42ti[t] ^ /^4t!ii+^[t\. Needless to say that 
all needed formulas have been also obtained in parallel for the Euclidean case. 



3 Applications to Higgs boson decay and Adler function 

First, we analyze the Higgs boson decay to a hb pair. Here we have for the decay width 

r(H^66) = -^Afffi?s(Afi) with RsiMjj) = mliMfi) RsiMfi) (13) 
4V27r 

and Rs{s) is the i?-ratio for the scalar correlator, see for details in [51[TU]. In the one-loop 
FAPT this generates the following non-power expansiorjfl: 

ns[L]=3mlA^i;^[L]+d! Y.^^tl\[L]y (14) 

I TI>1 ^ ) 

where ^(1) is the renormalization-group invariant of the one-loop m^(/i^) evolution 
ml{Q'^) = m^-^^ a'^°{Q^) with — 27o/^^o(5) = 1-04 and 70 is the quark-mass anomalous 
dimension (for a discussion — see in (llj). 

We take for the generating function P{t) the Lipatov-like model of [21 [H] with 
{c = 2.4, P = -0.52} 

Js_ n-i r(n-H)+/3r(n) ^ p,..m+£ ,/. 

It gives a very good prediction for with n — 2,3,4, calculated in the QCD PT \TU\ : 
7.50, 61.1, and 625 in comparison with 7.42, 62.3, and 620. Then we apply FAPT 
resummation technique to estimate how good is FAPT in approximating the whole sum 



^Appearance of denominators tt" in association with the coefficients dn is due to d„ normalization. 
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TZs[L] in the range L e [11, 13.8] which corresponds to the range Mh € [60, 170] GeV^ 
with Aq^d^ = 172 MeV and 2lf°''(m|) = 0.120. In this range we have {Lq = In(m2/A^)) 



7^s[L]=3m^l) ^2t^-[L] + ^((ai+ 



L + X5-—;5 
TrP5 ■ 



6-^1+1^0 



(16) 



with PJ^ {t) defined via Eqs. (llSp and (fT2|) . Now we analyze the accuracy of the truncated 
FAPT expressions 



N 



n=l 



(17) 



and compare them witli the total sum 7?.s[i] in Eq. (|16p using relative errors A^[L] = 
1 - ns[L;N]/ns[L]. In the left panel of Fig. H we show these errors for iV = 2, TV = 3, 
and iV = 4 in the analyzed range of L e [11, 13.8]. We see that already TZs[L;2] gives 
accuracy of the order of 2.5%, whereas 7?,s[L;3] of the order of 1%. That means that 
there is no need to calculate further corrections: at the level of accuracy of 1% it is quite 
enough to take into account only coefhcients up to d^. This conclusion is stable with 
respect to the variation of parameters of the model Ps {t) . 
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Figure 2: Left panel: The relative errors A^[L], TV = 2, 3 and 4, of the truncated 
FAPT, Eq. (I17|) . in comparison with the exact summation result, Eq. (I16|) . Right panel: 
Analogous relative errors A)^(Q^), N — 1, ... ,4, for the case of vector Adler function 
(solid line is for = 2, dashed — for A^ = 3, and long-dashed — for A^ = 4). 

Next, we analyze the Adler function of vector correlator. In the one-loop APT this 
generates the following non-power expansion: 

2?v[L] = l + <E§-^""W' (18) 

n>l 
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Here we use another Lipatov-like model for perturbative coefficients in the Nf — 4 region 

with c = 3.456 and (3 — 1.325, which gives a very good prediction for with n — 2,3, 4, 
calculated in the QCD PT [13]: 1.49, 2.60, and 27.5 in comparison with 1.52, 2.59, and 
27.4. To estimate how good is APT in approximating the whole sum [L] , we apply APT 
resummation approach in the range G [1.5,20] GeV^, corresponding to the Nj = 4 
value. Again we analyze the accuracy of the truncated APT expressions VviLiN] = 

^ + di J2n=i ^ ^n°^[-^ ^^'^ compare them with the total sum Pv[i]i obtained by 
resummation APT method, using relative errors A)^[L] = 1 — VviL; N]/'Dv[L]. In the 
right panel of Fig. [5] we show these errors for TV = 1, . . . , 4 in the analyzed range of Q^. 
We see that already V^IQ^; 2) gives the accuracy of the order of 0.05%, whereas taking 
into account higher-order corrections only worsen the accuracy: I?v((5^;3) provides the 
accuracy of the order of 0.1% and T>v{Q^; 4) — of the order of 0.2%. That means that 
the NLO approximation gives the best result and after that the series starts to reveal its 
asymptotic character. 

4 Conclusions 

In this report we described the resummation approach in the global versions of the one- 
loop APT and FAPT and argued that it produces finite answers in both Euclidean and 
Minkowski regions, provided the generating function P(t) of perturbative coefficients d„ 
is known. In the case of the Higgs boson decay an accuracy of the order of 1% is reached 
at N'^LO approximation, when term ^3^3 is taken into account, whereas for the Adler 
function D{Q^) we have an accuracy of the order of 0.1% already at N^LO (i.e., with 
taking into account d2A2 term). 

The main conclusion is: In order to achieve an accuracy of the order of 1% we do not 
need to calculate more than four loops and d^ coefficients are needed only to estimate 
corresponding generating functions P{t). 
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